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Abstract 

At hadron colliders invisible particles x can be inferred only through observation of the transverse component of 
the vectorial sum of their momenta - missing Et or MET - preventing reconstruction of the masses of their mother 
particles. Here we outline situations where prior prejudice about the event kinematics allows one to make the most 
of MET by decomposing it into its expected sum of transverse contributions, each of which may be promoted to a 
full four-momentum approximating the associated \- Such prejudice arises when all \ in the event are expected to 
be light and (anti-)parallel to a visible object, due to spin correlations, back-to-back decays or boosted decays. We 
focus on the last of these, with boosted semi-invisibly decaying neutralinos widely motivated in supersymmetry (in 
the presence of light gravitinos, singlinos, photini or pseudo-goldstini) , and demonstrate our simple method's ability 
to reconstruct sharp mass peaks from the MET decomposition. 



1 Introduction 

Missing transverse energy or MET is of great importance 
at hadron colliders: it is our only way of inferring the pres- 
ence of non-interacting (collider-) stable particles \- How- 
ever whenever two such particles are produced (which will 
always be the case if their stability is due to a Z 2 sym- 
metry) our observation only of the vectorial sum of their 
transverse momenta $ T = |j4 | = \p a ,T + Ph,r| thwarts 
reconstruction of masses in the decay cascad^J ending 
with 2x- Popular methods for searching for heavy par- 
ticles with partially invisible decays are transverse mass 
observables [l], Mti (2], razor analyses [3] and kinematic 
edges [4]. 

We start by asking under what circumstances can we 
explicitly access the missing momentum associated with 
each x particle separately. Clearly some feature of the 
rest of the event must suggest the correct decomposition 
of p into p Qj T + Pb,T- If there are two well-localised 
visible objects that we expect, from some prior prejudice 
about the kinematics, to be parallel or anti-parallel to the 
two unseen x particles, then we have two directions in the 
transverse plane to give us p a ,r and Pb,T- Furthermore 
we can add longitudinal components to each of these two 
transverse vectors to make them (anti-)parallel to their 
corresponding visible object in three dimensions, giving 
approximations for p Xa b . If x is much lighter than the 
particle produced in the hard scattering, i.e. at the start 
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1 2x could also be directly produced, giving a final state with, at 
leading order, no large transverse energy (visible or invisible). The 
universal possibility of hard initial state radiation allows essentially 
model-independent limits to be set on the direct production of new 
X particles from monojet and monophoton searches. Here we focus 
only on production of 1\ via a decay cascade. 



of the decay cascade, we can promote p Xa b to massless 
four-vectors; we will show that combined with the four- 
vectors for the visible decay products, a strong mass peak 
for the initial particles can be reconstructed. 

2 Motivation 

Parallel or anti-parallel visible and missing energy is not 
worth considering only for its ease: it can arise in many 
circumstances. Spin correlations may make x particles ap- 
proximately (anti-) parallel to other particles. Two-body 
decays of particles P not boosted in the lab frame, such as 
the majority of those produced directly, are usually back- 
to-back: therefore in 2P — > 2x+2vis, each x is often nearly 
anti-parallel to one of the 'vis'. (However this simple ex- 
ample lends itself well to numerical optimisation over all 
possible momenta for the 2x, i.e. to use of Mr 2 -) 

Of particular interest is when each x is produced to- 
gether with visible energy from the decay of a boosted par- 
ticle. This will arise whenever a) directly pair-produced 
particles are appreciably heavier than whatever they decay 
into in the first step of the cascade, and b) x are created 
following two or more steps. Together these points imply 
that each of the two 'sides' of the event (separated ac- 
cording to the mother particle) contains an intermediate 
particle which is boosted: the visible object(s) and x it 
ultimately decays to will be collimated. 

For some examples, consider the quintessential super- 
symmetry (susy) decay of a pair-produced squark to a 
hard jet and light neutralino : q — > q + N\ (wc denote 
the susy neutralinos by N± to avoid confusion with our 
generic invisible particle x). There are many reasons why 
we might expect N\ to be unstable, decaying to visible 
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energy and a lighter, neutral, collider-stable particle - the 
latter could be: 

• a gravitino G, if susy breaking is mediated at a low 
scale, i.e. some form of gauge mediation. A low me- 
diation scale is motivated by electroweak naturalness 
and an automatic solution of the susy flavour prob- 
lem. See [5] for a review and [6j for a comprehensive 
list of possible collider signatures. 

• a pseudo-goldstino G\ if more than one hidden sec- 
tor breaks susy, as may occur in string theory or 
quiver gauge theories [7j[8]- See [9] for the collider 
phenomenology. 

• a smg lincg S, if the Minimal Supersymmetric Stan- 
dard Model (MSSM) is extended with a gauge-singlet 
superfield, giving the Next-to-MSSM (NMSSM). This 
is motivated by the /Lt-problem of the MSSM, and also 
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by naturalness 10 . See 11 12 
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for the modified collider signals. 



• a new photino 7 , if the MSSM is extended with 
one or more extra U(l) gauge symmetries, as is com- 
monly expected to arise from string compactifica- 



tions. See 14 for a discussion of collider prospects. 



In the nomenclature of 
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N\ here is the Lightest Ordi- 
nary Supersymmetric Particle (LOSP). All of these ex- 
amples have some other particle as the true LSP, and 
so a charged or coloured susy particle could be lighter 
than Ni and take its place as the LOSP in the cascade 
q visi + (LOSP) -t visi + (vis 2 + LSP), giving different 
visible energy. 

3 The Analysis 

We elaborate on the strategy outlined in Section[T]in terms 
of a concrete example to allow clearer references to the 
particles involved in the signal: we consider the classic 
gauge-mediation decajj^2q -> 2q + 2(Ni) -> 2q + 2(G + 7). 
The lightest ncutralino is typically expected to be consid- 
erably lighter than the squarks in this scenario, as the 
phenomenon of gaugino screening in the simplest mod- 
els makes the gauginos much lighter than the scalars (see 
e.g. 16 ) and renormalisation-group evolution tends to 
drive squark masses up and the bino mass down. This 



2 In this case the decay is not really q — > q + N\ — > • ■ • but 
q — > q+N2 —> Ni + . . ., since new photini/singlinos actually mix with 
the MSSM neutralinos. If N2 is mostly 'MSSM-like' (any mixture 
of Higgsino, wino and bino), and N± is mostly singlino or a new 
photino, then direct decay of q to Ni is suppressed relative to the 
two-step decay. 

3 A similar final state may arise from Universal Extra Dimensions 
[15] , though semi-invisibly decaying Kaluza-Klein photons from KK 
quark/gluon decays are not generally expected to be boosted; this 
will be important for our analysis. 



simple observation gives a powerful handle on the signal, 
as yet unexploited: the gravitinos and photons are nor- 
mally collimated. It is exploited as follows. 

1. Uniquely decompose $ T into p a ,T + Pb,T which are 
defined to be parallel, in the transverse plane, to the 
two hardest isolated photons. 

2. Promote p a ,b-.T to three-vectors p a & by adding the 
longitudinal components required to make them par- 
allel to each of the photons in three dimensions. 

3. Promote p a .b to massless four-vectors p^ b = 
(|Pa,i>|i Pa,&), giving approximations for the two grav- 
itino four-vectors. Adding each of these to the four- 
vector of the collinear photon gives massless approx- 
imations for the two neutralino four-vectors, r> ~ 

4. If each neutralino Ni- a ^ can be paired with the 'cor- 
rect' jet in the event j a ,6, then taking the invariant 
mass of each pair reconstructs the mass of the initial 
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Steps 1-2 above reconstruct the three-momenta of the 
two neutralinos in the same way as is done for the two 
t in U — y 2t — V e ± /j,T$ T with the collinear approxima- 
tion of |17| (an approximation which we also found to be 
useful for jet substructure in [18]). There, the two r four- 
momenta are added together to get the mass of the single 
mother particle; here the four-momenta of the two neu- 
tralinos are separately added to those of other visible par- 
ticles in the event to get the masses of two mother particles 
- step 4. 

Step 4 needs a criterion for the correct way to pair each 
reconstructed neutralino with one of the jets in the event, 
since the mass they should reconstruct is unknown. The 
correct jet is considered to be the one most closely re- 
sembling the quark produced in the same q — > Ni + q 
decay. Keeping only the two hardest jets, there are two 
arrangements - two ways of pairing each neutralino with a 
different jet. More generally one can consider the N hard- 
est jets in the event, giving N(N — 1) arrangements to 
choose from. Each squark is generally produced nearly at 
rest, therefore the neutralino and jet into which it decays 
are likely to be back-to-back; the jet is also expected to 
be hard, with an energy roughly half the squark's mass. 
Therefore one criterion is to pair the two neutralinos -/Vi ;clj {, 
with jets j a and jj, so as to make maximally negative the 
sum of dot products between the three-momenta of each 
neutralino and its jet: 



criterion a : 



(p.y. P, 



maximal 



If the pair-produced squarks are mass degenerate, this can 
also be exploited: the two reconstructed masses should 
coincide. This gives the second possibility for finding the 
right jets: 



criterion (3 : 
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minimal 
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Each criterion suggests the correct jets, defining two re- 
constructed masses A/ 2 .„ . = (pfi ) 2 . The max- 

I (_,(_., LL , U y 1 J\j a b J a .b ' 

imisation/minimisation above is not differential but dis- 
crete - the quantity is calculated once for each of the 
N(N — 1) arrangements of jets with neutralinos and only 
the largest /smallest is kept. It thus takes negligible com- 
putational time (indeed N = 2 is optimal in our example) 
and could be incorporated into a trigger. These two cri- 
teria are not specific to neutralinos and jets: they are rel- 
evant for final states where two objects need to be paired 
correctly with two other objects, both being the decay 
products of pair-produced particles (the second criterion 
also requires mass degeneracy of the two mother particles) . 

We consider a simplified model with squarks of the first 
two generations, a bino-like neutralino, and a gravitino 
with masses mq — 1.2 TeV, rtiff = 100 GeV, = 
1 eV respectively; this squark mass is at the edge of 
the strongest current constraints [19 . We calculate a 
full spectrum for this simplified model (all other super- 
partner masses are set 2 TeV) with SoftSusy 3.3.4 20 
and decay widths with Herwig++ 2.6.1 |2l] . We then 
follow two routes to get to observable distributions. In 



the first, MadGraph 5 1.5.5 22 supplies the matrix el- 



ements for disquark production; the subsequent decays, 
extra radiation, showering and hadronisation are done 
by PYTHIA 6 23 ; fast detector simulation is then per- 
formed with PGS 4 24 . In the second, Herwig++ is 
used to generate the complete event; jets are defined with 
Fast Jet 3.0.3 25 , and the final state objects analysed 
in the RIVET 1.8.1 framework [26] . Our kinematical anal- 
ysis - steps 1-4 with criteria a and /3 above - is then ap- 



plied. Code for doing this can be found at 27 



Basic cuts needed for the analysis are as follows: 

• At least two jets, clustered using the anti-kt algo- 
rithm | 28] with size parameter 0.4. Jet candidates 
are required to have pt > 30 GeV and \r]\ < 4.5. 

• At least two isolated photons with pt > 10 GeV. 
On the MadGraph-PYTHIA route, PGS handles isola- 
tion. On the Herwig route, we consider a photon 
isolated when the sum of transverse energy in a cone 
AR < 0.4 around the photon is less than 5 GeV. 

• A minimum and maximum azimuthal angular sep- 
aration between the two hardest isolated photons 



e < A0 7l72 < 7r — e with e 



0.01, since photons 



which are exactly (anti-) parallel in the transverse 
plane do not allow $ T decomposition. 
The missing energy vector y> T should lie in between 
the two photons in the transverse plane (i.e. inside 
the smaller of the two sectors delimited by the two 
photon directions). This ensures that the event has 
fi T corresponding to the ansatz of both gravitinos 
being parallel to their photons. With this cut the 
kinematics are always in the 'trivial zero' of the Mt2 
observable (see 29 



Decomposition of $ T of course requires f T ^ 0; in prac- 
tice this is always satisfied. We do not cut on $ T - we 
analyse this particular signal not to optimise the asso- 
ciated cuts but simply as a demonstration of the mass 
reconstruction technique. In our present example almost 
all events have $ T > 100 GeV and so a large require- 
ment could be placed as in existing searches (likewise for 
the leading jet and photon px which are typically hard in 
the signal.) Note that with a requirement for hard pho- 
tons and $ T there is typically very small background for 
new physics [30] and the priority is an observable that in- 
creases the visibility of the signal alone, ideally through a 
resonance. 

Fig. [TJ shows the results of the analysis. 
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Figure 1: The squark mass (m,q = 1.2 TeV) in the process 
pp^rlq^r 2q + 2(Ni) -> 2q + 2(G + ^) reconstructed with 
missing energy decomposition and neutralino-jet pairing 
as described in the text. The masses of the lightest neu- 



tralino and gravitino are m 



100 GeV, m, 



1 cV; the 



centre of mass energy is 8 TeV. Panels on the left (right) 
show the mass of the squark calculated from the leading 
(sub- leading) photon in each event. Upper (lower) pan- 
els pair jets with reconstructed neutralinos using criterion 
a (/?). The blue dashed line shows events generated by 
MadGraph and PYTHIA, with fast detector simulation per- 
formed by PGS; the red solid line shows events generated 
by Herwig++. 



As our analysis makes use of hard jets arising from the 
decay of signal particles, it could in principle be affected 
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by the (higher order) production of additional jets in the 
hard scattering. To investigate this we simulated 2q and 
2q+ ljet production and combined these consistently into 
a single sample using the MLM matching procedure [31] . 
The reconstructed mass distributions are essentially iden- 
tical to those of simple 2q production shown in Fig. [T] 
which follows from the fact that our method is designed 
to find the two jets that look most like they have been 
produced by the decay of the squarks, and other jets are 
discarded. 

Criterion a can also reconstruct the masses of pair- 
produced non-degenerate particles. In Fig. [2] it is used 
to analyse the same signal as previously but now with 
one squark from the first two generations having mass 
1.1 TeV and the other seven having mass 1.4 TeV. This 
unequal splitting is chosen to have large cross-sections 
for the production of two squarks of different mass (four 
lighter squarks and four heavier would merely result in 
a dominant production of the lighter four alone); never- 
theless production of two squarks of the same mass still 
has non-zero cross-section. Thus the distribution of the 
larger (smaller) of the two masses calculated for each 
event peaks strongly at 1.4 TeV (1.1 TeV) and weakly 
at 1.1 TeV (1.4 TeV), with the weak peak resulting from 
pair-produced degenerate squarks. 
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Figure 2: As Fig. [l]but with one squark from the first 
two generations having mass 1.1 TeV and the other seven 
having mass 1.4 TeV. The solid red (dashed blue) line 
shows the smaller (larger) of the two masses reconstructed 
in each event. Only criterion a for jet-neutralino pairing 
is used. Events are generated with Herwig++. 



4 Discussion 

The final state of the example from the previous section 
has two jets and two pairs of roughly collinear photons 
and gravitinos. The jet could be replaced by any other 
visible particle - 'visi' - the photon too - l vis2' - and the 



gravitino by anything invisible, \: we show this general 
topology in Fig. [3] Provided there are two semi-invisible 
decays which are boosted (or forced into (anti-)parallel be- 
haviour by spin correlations) the same analysis presented 
here should in theory have some potential for mass recon- 
struction. Of course if visi and vis2 are objects less clean 
experimentally than light-flavour jets and photons, such 
as b quarks or even combinations of particles, the proce- 
dure will be more difficult in practice. Searches for mass 
peaks in the manner presented, considering various differ- 
ent particle types for visi,2, could discover expected or un- 
expected resonances. Below, we outline how the method 
might be adapted as the topology is distorted and gener- 
alised further. 




Figure 3: The topology we consider: pair-produced par- 
ticles each decay into a visible Standard Model particle 
visi and a much lighter particle, which is thus boosted; 
this decays semi-invisibly into vis2 and x- 



A Less Boosted Intermediate. Collinearity of x & n d vis2 
relies on their common mother particle being boosted; as it 
becomes less boosted they become less collinear. We show 
this effect, and the decreasing sharpness of the mass re- 
construction that results, in Fig. [4] for our previous gauge- 
mediation example, is increased from 100 to 400 GeV 
for constant rrig = 1.2 TeV. If N\ is made heavier still, e.g. 
mfj^/rtiq — > 1, the increasingly lethargic neutralino gives 
a less collimated photon-gravitino pair; indeed the two are 
increasingly back-to-back, and most events fail to meet the 
requirement that p T be in between the two photons. 

More Decays Of The Intermediate. If vis2 is several par- 
ticles instead of the single photon 7 we considered, e.g. a 
lepton pair from a boosted Ni — > I PNi decay, by con- 
struction they will be collimated and the sum of their four- 
momenta can be used in place of p!^ in the analysis. 

More Decays Before The Intermediate. If the directly 
pair-produced particles decay to a boosted intermediate 
and two visible particles rather than one - via two on-shell 
steps or a three-body decay - then each visi in Fig.[3]is re- 
placed by two particles which are not collinear. Criterion a 
is then not applicable but criterion /3 is, albeit with greater 
combinatorial ambiguity from the need to pair each recon- 
structed neutralino with two other visible objects. In this 
scenario the boosted intermediate is also less boosted from 
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sharing its energy with more particles, making its semi- 
invisible decay less collimated. Despite these difficulties 
the method is reasonably successful: Fig. [5] shows events 
for a simplified model with pair-produced gluinos of mass 
1.2 TeV decaying to qqNi (q now denoting a quark of any 
of the three generations) with the 100 GeV neutralino de- 
caying to 7G. Neutralino-jet pairing is performed with 
criterion j3 generalised in the obvious way to include four 
jets rather than two. 




0.4 0.8 1.2 1.6 2 0.4 0.8 1 .2 1 .6 2 

M rec [TeV] M rec [TeV] 




0.5 1 1.5 2 2.5 0.4 0.8 1.2 1.6 2 

AR(tG) M rec [TeV] 

Figure 4: As Fig. [T] holding the squark mass at 
1.2 TeV while varying the neutralino mass: = 
100, 200, 400 GeV are shown with red solid, green dashed, 
and blue dotted lines respectively. The greater , the 
less collinear its photon and gravitino daughters become, 
as shown by AR(jG) (averaged between the two 7G pairs) 
in the left panel. This worsens the mass reconstruction: 
the right panel shows one of the two masses found using 
one of the two jet-neutralino pairing criteria (all four quan- 
tities behave similarly - see Fig. [IJ . Events are generated 
with Herwig++. 



Other Combinatoric Complications. If visi = vis 2 , e.g. 
if in our former example photons were replaced by jets or 
jets by photons (but not both of these at once), then there 
would be a combinatoric ambiguity not just in pairing the 
reconstructed boosted intermediate with the correct visi 
but also in which two particles define the initial p de- 
composition directions. The requirement that p be in 
between the two visible particles onto which it is decom- 
posed eliminates some of the possible decomposition con- 
figurations; for the rest, criterion /3 can be generalised to 
be an optimisation over decomposition configurations as 
well as pairing possibilities. 

More Than Two Invisible Particles. With a third x m 
the final state which is expected to be (anti-)parallel to one 
of the first two, our ansatz for the topology still contains 
only two invisible directions and we can uniquely decom- 
pose the observed MET. If the two invisible particles that 



Figure 5: The gluino mass (rrig = 1.2 TeV) in the 
process pp -> 2g -> Aq + 2(Ni) -t Aq + 2(G + 7) 
reconstructed with missing energy decomposition and 
neutralino-jet pairing criterion /3 as described in the text. 
The masses of the lightest neutralino and gravitino are 
m ] y i — 100 GeV, = 1 eV (the masses of other parti- 
cles are set at 2 TeV); the centre of mass energy is 8 TeV. 
Panels on the left (right) show the mass of the gluino 
calculated from the leading (sub-leading) photon in each 
event. The blue dashed line shows events generated by 
MadGraph and PYTHIA, with fast detector simulation per- 
formed by PGS; the red solid line shows events generated 
by Herwig++. 



are (anti-) parallel have come from the decay of the same 
particle, we only need to know the sum of their momenta 
and so we can reconstruct the mass as before. However if 
they have come from the decay of two different particles, 
then we need their individual momenta for mass recon- 
struction; knowing only their sum, the masses we wish 
to calculate are underconstrained by one parameter. An- 
other possibility is 3x in the signal final state with three 
different expected directions: there are then three vectors 
in the transverse plane into which p can be decomposed, 
with any two of the three giving a unique decomposition. 
There are three ways to choose two vectors from the three. 
We may have the p T in between the two vectors in 0, 1 
or 2 of the three ways (neglecting the possibility of ex- 
act collinearity between p T and one of the vectors). If 
0, we veto. If 1, there is a unique decomposition. If 2, 
p T can be expressed as some amount of one of the de- 
compositions plus some amount of the other, with the two 
coefficients constrained to sum to unity: the masses we 
wish to calculate are under-constrained by one parameter. 
One response, not physically motivated, would be to veto. 
Which of these three cases (0, 1 or 2 of the possible de- 
compositions being acceptable) we have will vary on an 
event by event basis. 
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5 Conclusion 



References 



Partially invisible decays occur in many extensions of the 
Standard Model, being more or less omnipresent in mod- 
els with dark matter candidates. However when missing 
energy arises from two invisible particles, reconstruction 
of mass peaks - and hence discovery - is much more dif- 
ficult, even though Standard Model backgrounds may not 
be dominating the signal. We consider the case where 
there are two preferred directions for the two invisible 
particles, which may arise due to spin correlations, back- 
to-back decays or boosted decays. The observed missing 
energy vector can then be decomposed into components 
along each of those directions; these components approxi- 
mately describe the two invisible particles and may allow 
reconstruction of the mass of their mother particles. 

Pairs of collimated semi-invisible decays together with 
jets occur in supersymmetry when heavy squarks or 
gluinos decay to light MSSM-like neutralinos, which de- 
cay in turn to gravitinos, pseudo-goldstini, singlinos or 
photini. We considered as a concrete example the gauge- 
mediation style decay of 1.2 TeV squarks or gluinos to jets, 
photons and gravitinos. When the mass of the intermedi- 
ate neutralino is 100 GeV, the initial mass is reconstructed 
to 10% accuracy for roughly | of events passing the ba- 
sic cuts. Multiplying by the Prospino [32j[33] production 
cross-section and the acceptance - 20 fb x 0.5 for squarks, 
2.5 fb x 0.3 for gluinos - one would expect 0(100) events 
for squarks, O(10) for gluinos, in 30 fb -1 at 8 TeV. (These 
numbers of events inside the peak of course depend on the 
masses - when the mass of the neutralino approaches the 
mass of the squark/gluino, it is no longer boosted and a 
peak will not be reconstructed with this method.) 

We have discussed how the method is applicable to final 
states with particles different from those in the example, 
and outlined the limitations as it is applied to more gen- 
eral scenarios. The level of reconstruction is extremely 
encouraging and we hope that our results are an incen- 
tive for the experimental collaborations to investigate the 
feasibility of implementing such analyses. 
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